A lattice of double wells for manipulating pairs of cold atoms 
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We describe the design and implementation of a 2D optical lattice of double wells suitable for 
isolating and manipulating an array of individual pairs of atoms in an optical lattice. Atoms in 
the square lattice can be placed in a double well with any of their four nearest neighbors. The 
properties of the double well (the barrier height and relative energy offset of the paired sites) can 
be dynamically controlled. The topology of the lattice is phase stable against phase noise imparted 
by vibrational noise on mirrors. We demonstrate the dynamic control of the lattice by showing the 
coherent splitting of atoms from single wells into double wells and observing the resulting double-slit 
atom diffraction pattern. This lattice can be used to test controlled neutral atom motion among 
lattice sites and should allow for testing controlled two-qubit gates. 



o 
o 



> 
m 
o 

O 

o 



o 
o 



X 



Bose Einstein condensates (BEC) in optical lattices 
have proven to be an exciting and rich environment for 
studying many areas of physics, such as condensed mat- 
ter physics, atomic physics, and quantum information 
processing (see for instance IJ). Optical lattices are very 
versatile because they allow dynamic control of many im- 
portant experimental parameters. Dynamic control of 
the amplitude of the lattice has been widely used (e.g. 
0, 0, lji 0); recent experiments have used a state de- 
pendent lattice to dynamically control the geometry and 
transport of atoms in the lattice Q . Recently there have 
been several proposals for using optical lattices to per- 
form neutral atom quantum computation 0, IM 13 • With 
optical lattices it should be possible to load single atoms 
into individual lattice sites with high fidelity 10], and 
then to isolate and manipulate pairs of atoms confined by 
the lattice in order to perform 2-qubit gates. Loading of 
single atoms into lattice sites or traps was demonstrated 
by J:I, 12, 13, 14], but to date no neutral atom based 
trap can isolate and control interactions between indi- 
vidual pairs of atoms. While previous experiments have 
demonstrated the clustered en tang lement of many atoms 
confined by an optical lattice [l5j |. the unique ability to 
isolate and control interactions between pairs of atoms 
would allow for entanglement between just the pair of 
atoms. 

In this paper we report on a double well optical lat- 
tice designed to isolate and control pairs of atoms. The 
lattice is constructed from two 2D lattices with different 
spatial periods, resulting in a 2D lattice whose unit cell 
contains two sites. Within the pair, the barrier height 
and relative depths of the two sites are controllable. Fur- 
thermore, the orientation of the unit cell can be changed, 
allowing each lattice site to be paired with any one if its 
four nearest neighbors. The double well lattice is phase 
stable in that its topology is not sensitive to phase noise 
from motion of the mirrors. This lattice, in combination 
with an independent ID lattice in the third direction to 
provide 3D confinement, is ideal for testing many 2 qubit 
ideas, particularly quantum computation based on the 
concept of "marker atoms" @ and controlled collisions 
0. Among other applications, this lattice could be used 



for studying tunnel coupled pairs of ID systerns, inter- 
esting extensions to the Bose Hubbard model [laj, and 
quantum cellular automata [T^ . 

This paper is divided into six sections. In Section I we 
discuss the ideal structure of the lattice. Section II de- 
scribes several experimental issues which need to be con- 
sidered in order to experimentally realize an ideal double 
well lattice. Section III details the experimental real- 
ization of this lattice and a measurement of the impor- 
tant parameters. In Section IV we show the momentum 
components present in our lattice by mapping the lattice 
Bruillion zone. In section V we demonstrate the dynamic 
control of the properties and topology of the double well 
lattice by showing the coherent splitting of atoms from a 
single well into a double well. We summarize and present 
prospective applications in Section VI. 



I. IDEALIZED 2D DOUBLE- WELL LATTICE 

An ideal double-well lattice would allow for atoms in 
neighboring pairs of sites to be brought together into 
the same site, requiring topological control of the lattice 
structure. It has been shown |l8| that a D-dimensional 
optical lattice created with no more than D-l-1 inde- 
pendent light beams is topologically stable to arbitrary 
changes of the relative phases of the D-l-1 beams. This 
geometry is usually preferred since phase noise (e.g. that 
imparted by vibrational noise on mirrors) will merely 
cause a global translation of the interference pattern. To 
allow for topological control, a general double-well lat- 
tice will necessarily have more than D-l-1 beams, but it 
would be desirable to preserve the topological insensitiv- 
ity due to mirror-induced phase noise. To achieve vi- 
brational phase stability in a D-dimcnsional lattice made 
with more than D-l-1 beams, one can actively stabilize the 
relative time phase between standing waves Qj] 20]. Al- 
ternatively the lattice can be constructed from a folded, 
retroreflected standing wave, which forces the relative 
time phase between standing waves to be a constant 0| . 
Examples for a 2D case are shown in Fig. ^ 

In this paper we consider the latter design, shown in 
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FIG. 1: 2D lattices with four beams. (a) Lattices formed by 
interfering two independent standing waves must be actively 
stabilized to be topologically phase stable against phase noise 
caused by vibration of mirrors, (b) Lattices formed from a 
folded retroreflected beam have intrinsic topological phase 
stability. 



Fig. ^b. In this scheme, the same laser beam intersects 
the position of the atom cloud four times. The incoming 
beam with wave vector fci along x is reflected by mir- 
rors Ml and M2, and after traveling an effective distance 
di (where the effective distance includes possible phase 
shifts from the mirrors) returns to the cloud with wave 
vector k2. The beam is then retroreflected by M3, re- 
turning a third time with wave vector = — fe, having 
traveled an additional effective distance 2d2. Finally, it 
makes a fourth passage with — —ki, traveling again 
the distance di. The total electric field for this 2D 4- 
beam lattice is given by Re[i?(a;, ?/)e'"*], where 
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and r = XX + yy, 9 — kdi, (j) = fcc?2j k — I'k jX (A is the 
wavelength of the lattice light), and is the polarization 
vector of the i*'* beam. In the absence of polarization 
rotating elements and ignoring polarization dependent 
phase shifts from mirrors, 64 = e\ and 63 = 62. Since 
the beam retraces the same path, there are only two in- 
dependent relative phases between the four beams. As 
a result, the lattice is topologically stable to vibrational 
motion of Ml, M2, and M3; variations in d\ and d2 result 
in a simple translation of the interference pattern (21j . 

The potential seen by an atom in a field Re[i?e*'^*] is 
given by ?7 = — (l/4)£'* ■ a ■ E, where a is the atomic 
polarizability tensor 22] . In general, ct depends on the 
internal (angular momentum) state of the atom, having 
irreducible scalar, vector, and 2"'^ rank tensor contribu- 
tions with magnitudes a^, a„ and at, respectively. The 
scalar light shift, Us — — as|i?p/4, is state independent 
and directly proportional to the total intensity. The vec- 
tor hght shift, Uy = iay{E* x E) ■ F/i, depends on the 
projection of total angular momentum HE. It can be 
viewed as arising from an effective magnetic field whose 
magnitude and direction depend on the local ellipticity 



of the laser polarization, B^s ~ iav{E* x E). It vanishes 
for linearly polarized light. The total vector shift in the 
presence of a static magnetic field B is determined from 
the energy of an atom in the vector sum field B^s + B. 
The 2"'^-rank tensor contribution is negligible for ground 
state alkali atoms far detuned with respect to hyperfine 
splittings |22| . and we will ignore it in this paper. 

Consider the ideal situation with four beams of equal 
intensities {Ei — E) which intersect orthogonally {ki ■ 
k2 = 0). As a first case consider ei = y, 62 = x, where all 
the light polarizations are in-the-plane. We will refer to 
this configuration as the "in-plane" lattice. The spatial 
dependence of the electric field is given by the real part 
of 
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ly{x,y) = £:(^e'''^^ + e^(2«-y+20=.«-fe^))y 



where 9xy and 



'xy are the path length differences for 
in-plane light taking into account that the path length 
difference could be polarization dependent. This gives a 
normalized total intensity of 

Ixy{x,y)/h = 2cos{2kx - 29xy - 2(j)xy) 

+2cos{2ky + 2(j)xy) + 4: (2) 

where Iq is the intensity of a single beam. Due to the 
orthogonal intersection {ki ■ k2 — 0, etc.) and the or- 
thogonality of the polarizations between ki and k2 etc., 
the resulting four beam lattice is the sum of two inde- 
pendent ID lattices. As shown in Fig. this creates a 
2D square lattice with anti-nodes (and nodes) spaced by 
A/2 along x and along y. Since the four beam intensities 
are equal, the lattice forms a perfect standing wave, and 
the polarization is everywhere linear, although the local 
axis of linear polarization changes throughout the lattice. 
In this case the vector light shift vanishes, and the light 
shift is strictly scalar U{x,y) — —aseo\E{x,y)\'^/4:. Note 
from Eq. |2]that varying 9xy changes the relative position 
of the lattice formed by fci and k^, moving the lattice 
along X. The phase (pxy affects both ID lattices, shifting 
the combined 2D lattice along {x — y)/\/2. 

As a second case consider (ei — £2 — z), where all 
the light polarizations are out-of-the-plane. We will refer 
to this configuration as the "out-of-plane" lattice. The 
electric field is given by the real part of 
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where 9z and (j)z are the path length differences for out- 
of-plane light. In this case the intensity is not simply a 
sum of independent functions of x and y, but rather given 
by 

Iz{x,y)/lQ = 4:[coa {kx - 9z - (j)z) + cos {ky + (j)z)]'^ 
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FIG. 2: Calculated intensities for in-plane lattice (a) and the 
out-of-plane lattice (b). Cross sections taken on the white 
dashed line are shown below their respective plot; a cross is 
used to denote the origin in each plot. The in-plane lattice has 
the familiar cos^ profile typical of A/2 lattices, while the out- 
of-plane lattice has a cos^* profile and periodicity of A. The 
flat portion of the (b) cross section shows the intersection of 
two nodal lines. 
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As shown in Fig [2b, the added interference creates com- 
ponents at k in addition to the components at 2k re- 
sulting in a lattice spacing along x and y of A rather 
than A/2 (the lattice period along i -I- y is \/V2). In 
addition, the nodal structure changes in that there are 
nodal lines along the diagonals. In particular, every other 
anti-node of the in-plane lattice is at the intersection of 
two nodal lines in the out-of-plane lattice. The polariza- 
tion is everywhere linear along z, giving rise to a strictly 
scalar light shift. As with the in-plane lattice, varying 6z 
translates the out-of-plane lattice along x, and varying 
02 translates the lattice along [x — y)/^/2. 

A double well lattice is realized by combining the in- 
plane and out-of-plane polarizations. Since the polariza- 
tions of the two lattices are orthogonal, the total intensity 
is hot — I'xy + Iz, and the scalar part of the light shift 
is simply a sum of the light shifts from the in-plane and 
out-of-plane lattices. Electro-optic elements in the beam 
paths di and d2 can produce different phase shifts for 
different input polarization, allowing for control of the 
relative phases 59 = 9z — Oxy and 5(j) = (pz — 4>xy, while 
maintaining vibrational phase stability of the combined 
lattice. This combined lattice can have a vector light 
shift, since relative phase shifts between the two polar- 
izations allow for non-zero ellipticity, i{E* x E) 0. If 
both lattices are everywhere in time-phase {SO = or tt 
and S(t) — or tt), the vector shift vanishes. Otherwise, 
there is a non-zero, position dependent Bcff{x,y) which 
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FIG. 3: Adjustment of the phases 59 and S(j} allow for nearest 
neighbor pairing with all four nearest neighbors. marks 
the location of a lattice site located at the origin which can 
be paired with any of its four nearest neighbors (shown with 
O) depending upon the choice of phase: a) 59 = n/2,5(f> = 
-Tv/2 b) 59 = -tt/2,5(I) = 7r/2 c) 59 = -tv/2, 5(t) = d) 
59 = n/2, 5(t> = 0. 



lies in the x-y plane. 

Control of the phase shifts, d(j) and 69, and the rel- 
ative intensity, Ixy/Iz, provides the flexibility to adjust 
the double- well parameters: the orientation (which wells 
are paired), the barrier height, and the tilt. For in- 
stance, double-well potentials along the ai-direction can 
be formed by setting 5(f) = Q and 59 = 7r/2. Fig. O 
demonstrates how a site can be paired with any one of 
its four nearest neighbors. Control of the barrier height 
and of the tilt are shown in Fig. 0] 



II. REALISTIC 2D DOUBLE WELL LATTICE 

In the previous section we considered idealized lattices, 
making assumptions about the amplitudes, wave-vectors 
and polarizations of the beams in the lattice. In this sec- 
tion we discuss considerations needed to experimentally 
realize the lattices described above. 



A. In-plane lattice 

For certain applications, such as the realization of the 
Mott-insulator state [sll , we need a nearly perfect in-plane 
lattice, namely a square 2D lattice with little or no energy 
offsets between neighboring sites. There are three pri- 
mary sources of imperfections that affect the performance 
of the in-plane lattice: imperfect control of the input po- 
larization [ii ■ z = sin/3 ^ 0), imperfect alignment caus- 
ing the beams to be nonorthogonal, (fci • ^2 = sine ^ 0), 
and imperfect intensity balance among all four beams 

When trying to make a perfect in-plane lattice, if the 
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a) 1^/1= =9 b) 59=?t/3; 





FIG. 4: Cross sections of example double well potentials. 
Solid line represents the double well potential; dotted line 
shows the placement and amplitude of the out-of-plane lattice, 
a) The barrier height, labeled above by the quantity 7, of the 
double well can be adjusted by placing the out-of-plane lattice 
"in the barrier" and adjusting the ratio of Ixy/Iz- b) The 
"tilt" of the double well (the relative offset between adjacent 
sites) can be changed by adjusting 59 and S(j>. 



input polarization is tilted by an angle (3 with respect to 
the xy-plane, then there is a z component to the light. 
The result is a contamination of the in-plane lattice by 
an out-of-plane lattice that modulates the lattice depth 
with a periodicity of A (Fig. [Sji). Neighboring sites will 
experience an energy shift /SXJ — 4L/o sin^ (/3) where C/q 
is the depth of a /? = in-plane lattice. Since AC/ scales 
as 0^ for small /3, the in-plane lattice is fairly tolerant to 
small rotations of the input polarization. For example, a 
misalignment of 10 mrad will cause a 0.04% modulation 
of the trap depth. 

The more stringent demand for minimizing site-to-site 
offsets of the in-plane lattice is the orthogonality of the 
two standing waves. If fci • fc2 = sine 7^ 0, standing 
waves ki , ^4 and k2 , fcs have nonorthogonal polarization 
and give rise to an interference term in the total intensity, 
thus causing an energy offset between neighboring sites 
given by MJ ^ 4C/oe for small e (Fig. [SJi). This imperfec- 
tion has the same effect as imperfect input polarization, 
but is harder to minimize since it scales linearly with e. 
For example, a misalignment of 10 mrad will cause a 4% 
modulation of the trap depth. We describe below how to 
control both imperfections. 

The third source of imperfections for the in-plane lat- 




FIG. 5: Lattice imperfections causing a) modulation of the 
lattice depth, Af/, between neighboring sites and b) state de- 
pendent modulation of the barrier height by a polarization 
lattice. In b) the solid line is the cross section of the inten- 
sity lattice; the dashed line is the cross section of the state 
dependent lattice resulting from unbalanced beam intensities. 
Atoms in the ground state of each well are shown schemati- 
cally. 



tice is the intensity imbalance between the four beams. 
Experimentally, intensity imbalance can arise from re- 
flection and transmission losses along the beam path as 
well as from unequal beam waists at the intersection 
psf. In general, light imbalance breaks the symmetry 
between the x and y direction, which removes the degen- 
eracy between the vibrational excitations along x and 
y. Typically, this does not adversely affect the lattice. 
We also note that since the beam experiences the same 
losses while traversing di each time, then for equal beam 
waists EiEi = E2E3, and the losses do not produce well 
asymmetries. 

A more important consequence of intensity imbalance 
is that the total field is not everywhere linearly polarized, 
but rather has some ellipticity, 

^{E* X E) = [iEiE2 - EsEi) sm{kx + ky - 6) 
+ {E1E3 - E2E4) sm{kx -ky-e~ 2<j))] z. (4) 

This causes a state dependent spatially varying vector 
light shift with period A, even in the absence of the out- 
of-plane lattice. As evident from Eq. 0] for perfect inten- 
sity balance the ellipticity will vanish, resulting in purely 
linear polarization. Comparing Eq. [21with Eq. 2|onc can 
see that the phase of the polarization lattice is spatially 
out of phase with the intensity lattice (see Fig. re- 
sulting in a state dependent barrier height between lattice 
sites, with relatively little modification of the potential 
near the minima. 



5 



B. Out-of-plane lattice 

In general, the structure of the out-of-plane lattice is 
fairly robust against the three imperfections mentioned 
above. A minor consequence of field imbalance is the 
possible disappearance of perfect nodal lines. One finds, 
for example, that at the position of the nodal line inter- 
section, the intensity becomes 

Iz.mm — [El — E2 — E3+ Eif (5) 

where eg is the electric constant (permittivity of free 
space), and c is the speed of light in vacuum. There are 
exact nodes when Ei+ E4 = E2 + E3, and this condition 
is trivially satisfied when the light fields are balanced. 
As with the in-plane lattice, degeneracies of vibrational 
excitations are lifted when the intensities are imbalanced. 



C. The double- well lattice 




A composite in-plane and out-of-plane lattice can be 
made by adjusting the angle f3 to control the admixture 
of the two components. For the combined lattice, the 
consequences and control of imperfections are similar to 
the in-plane lattice. With the added flexibility to control 
intensity and relative phase, we can in fact use (3 and S9 
to compensate for e 7^ (at least for a given magnetic 
sub- level). The vector light shift for an intensity imbal- 
anced double-well lattice is somewhat more complicated 
(yet easily calculable), having position dependent ellip- 
ticity along x, y and z. Many experiments are carried 
out in the presence of a spatially uniform bias field -B, 
so that the total field seen by the atoms is given by the 
vector sum Stot — B + Bcs- For \B\ >> |-Boff |, the di- 
rection of the quantization axis remains nearly constant 
along B throughout the lattice. The magnitude of the 
state dependent shift in this limit is proportional to 




and only the component of B^s along B contributes 
to the potential. The ability to adjust the direction 
of B provides significant flexibility in designing state- 
dependent potentials, and allows for state dependent mo- 
tion of atoms between the two sites of the double-well. 



III. IMPLEMENTATION 

This double well lattice was implemented on an appa- 
ratus described elsewhere 0. ^^Rb Bose Einstein con- 
densates are produced in an ultra high vacuum glass cell. 



FIG. 6: Schematic of the experimental implementation of 
the 2D double well lattice made from a single folded, retrore- 
flected beam. Mirrors Ml and M2, lenses LI and L2, and 
E0M6' are mounted on a fixed plate. 



We use RF evaporation to make BECs with « 200, ODD 
atoms in the -F = 1, mp = — 1 hyperfine state. The 
BEG is confined in a cylindrically symmetric magneto- 
static trap with ix}±/2TT = 24 Hz and ijJ\\/2tt = 8Hz. 
The Thomas- Fermi radii of condensates are « 15/xm and 
« 40//m respectively, with mean-field atom-atom inter- 
action energy approximately 500 Hz. Atoms in the BEG 
are then directly loaded into the "tubes" created by the 
2D double well lattice potential. The lattice beams are 
derived from a continuous wave (GW) TiiSapphire laser 
with A = 810 nm, detuned far from the Dl (795 nm) 
and D2 (780 nm) transitions in ®^Rb. On average 2600 
in-plane lattice sites or 1300 out-of-plane lattice sites 
(tubes) are filled with approximately 80 and 160 atoms 
per site respectively. Due to the tight confinement, the 
mean-field energy is much larger in the tubes than in the 
magnetic trap, as much as 7 kHz. During our experi- 
ments the magnetic confining potential is left on. 

The experimental schematic of the double well lattice 
is shown in Fig. |H| An acousto-optical modulator (AOM) 
provides rapid intensity control of the lattice light. The 
lattice light is coupled into a polarization maintaining 
fiber to provide a clean TEMqo spatial mode. A Glan- 
Thompson polarizer after the fiber creates a well defined 
polarization in the xy-plane. The light is folded by plane 
mirrors Ml and M2 then retroreflected by concave mirror 
M3. Lenses Lq, Li, and L2, in the input beam and after 
Ml, M2 respectively provide a weak focus (all four beams 
have 1/e^ beam radius of w 170 /im) at the intersection 
of the four beams. A 1 cm thick optical flat after L2 is 
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FIG. 7: Experimental images of atom diffraction from a 3 /is 
pulse of a) the in-plane lattice lattice and b) the out-of-plane 
lattice after 13 ms time-of- flight. 



used to translate the beam with wave vector ^2 without 
changing the angle of ^2 relative to ki . Mirror M3 images 
the intersection point back onto itself. 

Three electro-optic modulators (EOMs): EOM/3, 
E0M6', and EOM(/), control the topology of the lattice. 
EOM/3 is aligned with its fast axis orientated 45° rela- 
tive to the axis of the Glan-Thompson polarizer, allow- 
ing for control of the angle f3, which determines the ratio 
Ixy/Iz = cot2/3. EOM^ and EOMcj) are ahgned with 
their fast axes in the xy-plane allowing for control of the 
differential phases 66 and Scj) respectively. For these ini- 
tial experiments EOM0 was not implemented. 

Li, L2, Ml, M2, and E0M6' are located on a fixed 
plate. A preliminary alignment of the optics on the fixed 
plate was performed before installation on the BEG ap- 
paratus. In particular. Ml and M2 were first aligned 
using a penta-prism, and then lenses Li and L2 were in- 
serted and aligned to minimize deflections. The entire 
plate was mounted next to the BEG apparatus and the 
input lattice beam, fci, was aligned to pass through the 
center of Li and L2. With this technique we measured 
that we were able to initially align the beams so that the 
intersection angle deviated from orthogonality by only 
I e 1=7 mrad. 

Galibration of the in-plane lattice depth is achieved 
by pulsing the lattice and observing the resulting mo- 
mentum distribution in time-of- flight (TOF) [2Jj. This 
atomic diffraction pattern reveals the reciprocal lattice 
of the optical lattice. Diffraction from the perfect in- 
plane lattice has momentum components at multiples 
of zL2hkx and ±2hky, while diffraction from the out-of- 
plane lattice has additional components at multiples of 
±\/2hk{x ± y) . The diffraction patterns for both lattices 
after 13 ms TOF are shown in Fig. |7| For 120 mW and 
at A = 810 nm, we measure an average lattice depth of 
Uo = 40Efl (Er = h^k'^/{2m) = h x 3.5 kHz, m is the 
Rubidium mass) in each of the independent ID lattices 
making up the in-plane lattice. As seen in Figl^h, we cal- 
culate that the out-of-plane lattice is four times deeper 
than the in-plane lattice for equal intensity. 

Pulsing the lattice is a useful method for determining 
the average in-plane lattice depth, but this method dis- 



closes little information about variations in depth AU 
between adjacent sites of the in-plane lattice (such as 
variations caused by /3 7^ and/or e 7^ 0). On the 
other hand, the ground state wave function of the in- 
plane lattice is sensitive to AC/, and we can use this to 
make /3, e « 0. Information about the ground state can 
be revealed by adiabatically loading the atoms into the 
ground band of the lattice 0, quickly switching off the 
lattice, and observing the atomic momentum distribution 
in TOF. In this technique the lattice must be turned on 
slowly enough to avoid vibrational excitation but quickly 
enough to maintain phase coherence among sites; for our 
parameters the timescale for loading is ~ 500 /xs. (Note 
that band adiabaticity is more complicated when we com- 
bine the in-plane and out-of-plane lattices to create a 
double well lattice since the tunnel couplings and tilt be- 
tween double well sites can create situations where band 
spacings are very small.) For a small but nonzero AU, 
this timescale is not adiabatic with respect to tunnel- 
ing between neighboring sites. In this way atoms are 
loaded into every site, even though the true single parti- 
cle ground state fills every other site. Therefore, atoms 
are not in an eigenstate of the potential, and the atomic 
wavefunction evolves in time. In such a lattice poten- 
tial, pictured in Fig. (Sji, atoms in adjacent sites acquire 
a differential phase, AUt/h. The ground band diffrac- 
tion pattern changes in time as the atoms are held in the 
lattice and the differential phase is allowed to "wind up" . 

To quantify the "ground band diffraction" patterns, we 
define a variable G given by 



G 



Nik - N2k 
Nik + N2k 



(7) 



where is the number of atoms with momentum com- 
ponents ±2hkx and ±2hky, and Nik is the number of 
atoms with momentum components ±\/2hk{x ± y) (see 
Fig. IHt). G is normalized so that the value G — —1 
corresponds to a diffraction pattern containing only mo- 
mentum components associated with the in-plane lattice. 

We use the ground band diffraction to set the input 
polarization to /3 = by observing the dependence of the 
diffraction pattern on the differential phase shift 69 at a 
fixed time. For (3 = the light has no out-of-plane com- 
ponent so that changing 66 with E0M6' does not change 
the topology but merely translates the lattice. The cal- 
ibration of EOM/3 is done by finding the condition in 
EOM/3 which ehminates the effect of EOM0, this corre- 
sponds to /3 w 0. In practice for a setting of EOM/3, 
several ground band diffraction images are analyzed at 
different values of 66. EOM/3 is then adjusted until scans 
of 66 produce no noticeable difference in the diffraction 
pattern. 

Sample data for the calibration of /3 is shown in Fig. 
This method for determining /3 = is convenient 
because it is independent of other lattice imperfections, 
in particular this method does not rely on e = 0. For 
example the optimal /3 for the data shown in Fig. 
occurs for GwOt^— 1. G«0 has no experimental 
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FIG. 8: a) Schematic of the momentum components that con- 
tribute to G = {Nik - N2k)/{Nik + N2k)- N2k is the sum of 
atoms in the momentum components designated with filled 
circles, and A'^ife is the sum of atoms in momentum como- 
nents designated with open circles, b) Calibration of EOM/3: 
/3 ~ 52 mrad (open circles), P ~ 34 mrad (triangles), and /3 « 
mrad (filled circles). 



significance; it depends only on the time that the atoms 
were held in the lattice. A perfect in-plane lattice would 
have G = — 1 for all values of S9 at all values of time . 
The fact that G 7^ — 1 for /3 = indicates the presence of 
momentum components at ±^/2hk{x ± y) due to e ^ 0. 
With this method we can set /3 = to zero within 17 
mrad, placing an upper limit on AU/Uq ~ 0.1%. 

After setting f3 — we determine e by looking at the 
time dependence of the ground band diffraction pattern. 
We adiabatically load the lattice in the method described 
above, then we observe the time oscillations in the ground 
band diffraction pattern varying between a diffraction 
pattern with G = — ltoG = -|-l. From the time evolu- 
tion of G (see Fig. O, we extract the misalignment of the 
intersection angle, e = AU/AUq. The data (open circles) 
in Fig. iniwere fit to an exponentially decaying sinusoid 
(solid line). 

It is interesting to note the substantial decay in the 
amplitude of the oscillations in G shown in Fig. Et^, and 
the reduced rate of decay in Fig. We do not fully 
understand this damping, or the reason why the damp- 
ing is much less for the improved e. Inhomogeneities 
in the lattice depth due to the Gaussian nature of the 
lattice beams are not large enough to account for the 
decay. However, factors such as mean field effects, tun- 
neling, and misalignments between the lattice beams and 
the magnetostatic trap could contribute to the damping. 
Regardless of the cause of the decay, we can use this 
method and the data shown in Fig. to calculate and 
improve e. 

From the fit to the time evolution of G we extract an 
oscillation frequency, which can be used to calculate e. 
We calculate |e| after the initial penta-prism alignment 
to be 7 mrad ± 0.2 mrad (Fig. |5^); the energy differ- 
ence between neighboring sites of a 40 lattice was 3.9 
kHz ± 100 Hz. We reduced e by adjusting M2 and the 
optical flat in order to change the angle of k2 while keep- 
ing the beam aligned on the BEG, then remeasured the 
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FIG. 9: Time dependence of the value G characterizing the 
diffraction patterns for atoms loaded into lattices with a small 
offset energy AU caused by e 7^ 0. Open circles are data 
points; solid lines are a fit to the data using a exponentially 
decaying sinusoid. The frequency of the oscillations given 
from the fit is inset in each image. From this frequency we 
determine e: a.) 1/ = 3900 Hz corresponds to e ~ 7 mrad, and 
h) 1/ — 775 Hz corresponds to e ~ 1.4 mrad. The data in a) 
was taken after the initial penta-prism alignment; the data 
in b) was taken after several iterations of measuring the fre- 
quency and then realigning the beams to further improve the 
angle. Schematics of the diffraction patterns corresponding 
to different values of G at different times are shown in the 
insets. The initial phase of G in a) and b) is arbitrary; it 
depends only on how much phase has been wound up during 
the loading time. 



oscillation frequency of G. After several iterations of re- 
alignment and measurement, we improved the alignment 
to |e| =1.4 mrad ± 0.2 mrad, which corresponds to an 
energy offset of 775 Hz ± 70 Hz for a 40 E^ lattice (Fig. 
Et). For a 10 E^ lattice the energy offset would be < 200 
Hz. It is clear given the signal-to-noise ratio in Fig. [SJd 
that, if required, the angle could be further improved. 

We estimate the amount of polarization lattice from 
the measured intensity imbalance of the four beams. The 
losses are due to imperfect anti-reflective coatings on 
optical elements and the uncoated glass cell. The rel- 
ative depth of the polarization lattice is a function of 

ay\E* X E\/ ^ttsl-Bp^ . For far-off- resonant traps the 
ratio {ay /as) becomes small for the 5si/2 ground state 
of ^''Rb H^ll^l, thus decreasing the polarization lattice 
depth. From Eq. (HJ the size of the vector potential Uy 
depends on the size and orientation of the bias field. For 
our measured intensities, I2 = 0.85/i,/3 — 0.81/i,and 
I4 = 0.70/1 with A = 810 nm, we estimate for the purely 




FIG. 10: Experimental images after 13 ms TOF of atoms 
filling the first Brillouin zone for a) the in-plane lattice and 
b) the out-of-plane lattice. The shapes of the BZs reflect the 
momentum components in each lattice. 



in-plane lattice that the maximum ratio U^/Ug in the 
barrier is w 8%. In a combined in-plane and out-of-plane 
lattice, the vector shift can lead to a state dependent tilt. 
For present experiments the vector shifts are not impor- 
tant, but in future experiments this could be useful to 
produce state dependent tunnel couplings and state de- 
pendent motion. 
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FIG. 11: a) Single-slit diffraction pattern resulting from a 
loss of phase coherence among out-of-plane lattice sites, b) 
Double-slit interference pattern caused by coherence between 
atoms within a particular double well but not among the en- 
semble of double wells. 



IV. VISUALIZING THE BRILLOUIN ZONE 

After minimizing the imperfections in the lattices, we 
can look at the Brillouin zones (BZ) for each of the two 
lattices. We load atoms into the lattice in 100 ms, a 
timescale that is slow with respect to both vibrational 
excitations and atom-atom interaction energies so that 
atoms homogeneously fill the lowest band. The lattice 
is then turned off in 500 /zs, mapping the atoms' quasi- 
momentum onto free particle momentum states 0, 0, 
I25I . Atoms that occupied the lowest energy band 
of a lattice will have momentum contained in the first 
BZ of that lattice. The mapped zones for both the in- 
plane and out-of-plane lattice arc shown in Fig. ^| As 
expected the bands are different for the different lattices. 
This is clear evidence that we have two distinct lattices 
with distinct momentum components. 



lattice. The width of the single-slit pattern is inversely 
proportional to the Gaussian width of the ground state 
wave function in each lattice site @. 

To demonstrate the coherent splitting of atoms, we 
start with the ground-loaded out-of-plane lattice, then 
dynamically raise the barrier to transfer the atoms into 
the symmetric double well lattice. The barrier is raised 
in 200 /iS by increasing the ratio of Ixy/Iz with EOM/?, 
while EOM6I is set to 56 7r/2 '30'|. This timescale is 
chosen to be slow enough to avoid vibrational excitations 
but fast enough to maintain phase coherence within a 
double well. Since there is no phase coherence from one 
double well to another, the resulting momentum distribu- 
tion is an incoherent sum of essentially identical double- 
slit diffraction patterns (shown in Fig. Wlh ) from each of 
the wavefunctions localized in individual double wells. 



VI. CONCLUSIONS 



V. DYNAMIC CONTROL OF THE DOUBLE 
WELL LATTICE 

As an example of the dynamic control of the double 
well lattice, we demonstrate coherent splitting of atoms 
from single wells into double wells. Initially, we load into 
the ground band of the out-of-plane lattice. The time 
scale for loading (100 ms) is sufficiently slow to ensure 
dephasing of atoms in neighboring sites. If at this point 
in time we suddenly turn off the lattice and allow 13 ms 
TOF, we observe a single, broad momentum distribution, 
shown in Fig. Illb . Since atoms on separate sites have 
random relative phases, this distribution is an incoherent 
sum of "single-slit" diffraction patterns from each of the 
localized ground state wavefunctions in the out-of-plane 



The ability to isolate individual atoms in controllable 
double well potentials is essential for testing a variety 
of neutral atom based quantum gate proposals. Two- 
qubit gate ideas typically involve state dependent motion 
0, II] or controlled state dependent interaction , but 
nearly all require the ability to move atoms into very 
near proximity (e.g. into the same site) and subsequently 
to separate them. The flexibility and dynamic control 
of the double well lattice can be used to demonstrate 
and test motion of atoms between wells. Furthermore, 
state dependence of the barrier height can be used for 
state dependent motion between wells, allowing for the 
possibility of 2-atom gates. 

In conclusion we have demonstrated a dynamically 
controllable double-well lattice. The geometry of this 
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lattice is topologically phase stable against vibrational 
noise, yet allows topological control of the lattice struc- 
ture. The design of the double well lattice allows for 
flexible real-time control of its properties: the tilt and 
the tunnel barrier between sites within the double well. 
In addition, the orientation of the double well can be ad- 
justed so that a site can be paired with any one of its four 
nearest neighbors. We have described technical issues 
and imperfections of the double well lattice, and we have 
presented techniques to minimize the imperfections. We 
have demonstrated dynamic control of the double well 
lattice by showing the coherent transfer of atoms from 
single wells to double wells. In the future, the double 
well lattice presented here could be used for applications 
in quantum computation and quantum information pro- 



cessing, as well as studying interesting extensions of the 
Bose-Hubbard model, such as the emergence of super- 
solids and density waves ,16] . 
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